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METEOROLOGICAL ASPECTS OF STORM SURGE GENERATION@ 


D. Lee Harris! 
(Proc. Paper 1859) 


INTRODUCTION 


the coastal floods, produced by Hurricane AUDREY in June 1957 claimed 
e than 300(1) lives and destroyed more than 100 homes. Even in this en- 
tened year of 1957, when hurricane warnings are broadcast by radio and 
vision long before the storm arrives, the high tides produced by the storm 
inue to claim more lives than all other aspects of hurricane behavior. Al- 
igh many of the newcomers in the Cameron area evacuated the threatened 
on soon after the hurricane warnings were announced, many of the old- 
“rs remained because they had weathered previous hurricanes with little 
culty and they did not expect this one to be more severe than the others. 
w facts concerning hurricanes and their effects on sea level, which if fully 
srstood, may prevent future disasters of this kind, are presented in this 
r. 

‘his paper is concerned primarily with the meteorological aspects of 

m surge generation by hurricanes. The physical aspects of hurricane 

m surge behavior along the Atlantic and Gulf coasts of the United States 

» been discussed in greater detail in Harris.(2 Zetler 3) gave an extensive 
tment of the effects of hurricanes on sea level at Charleston, S. C. 

field and Miller(4) considered the effects of several recent Atlantic Coast 
‘icanes. Hubert and Clark(5) published charts showing the high water 

ks established after 16 Atlantic and Gulf coast hurricanes. Reid(6) de- 
ped a storm tide index based on the size and intensity of the hurricane and 
slope of the continental shelf. 


; Discussion open until May 1, 1959. To extend the closing date one month, a 
ritten request must be filed with the Executive Secretary, ASCE. Paper 1859 is 
art of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
merican Society of Civil Engineers, Vol. 84, No. HY 7, December, 1958. 
resented at the ASCE Hydraulics Conference in Cambridge, Mass., 
ugust, 1957. 

fice of Meteorological Research, Storm Surge Research Unit, U. S. 
Jeather Bureau, Washington, D. C. 
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Hurricane Models 


The wind and pressure distributions of a hurricane are never known wi 
desired accuracy before the storm strikes land. Hence, in forecasting stor 
surges and in the design of engineering works to be constructed as prot 
from some future hurricane, it is necessary to use a storm model. Such 
models can be constructed in a variety of ways depending on the use for wh 
they are intended. For storm surge work it is desired to describe as ac 
rately as possible the distribution of atmospheric pressure and sustzined 
speeds near the surface of the sea. 

The pressure in many hurricanes can be represented with a close 
mation by an equation of the type 


Foe Pree Pew. (r/R), 


where P,, = the pressure just outside the storm, 
Po = the pressure at the center of the storm, 
R =some characteristic radius of the storm, 
r = distance from the center of the storm, 
F(r/R) indicates a function of r/R. 


The radius of the eye, the distance from the center to the zone of. maxim 
winds, and the distance from the center to the outermost closed isobar hav 
all been suggested for use as the characteristic radius. The proper choice 
probably depends on whether it is desired to obtain the maximum accuracy 
the high wind speed zone or for the entire storm. (Py, - Po) is intended as 
measure of the intensity of the storm, but P, is likewise difficult to define 
explicitly. Some workers prefer to use the pressure of the outermost clo: 
isobar, still others prefer some other definition. Virtually all workers in 
this field use the best obtainable estimate of the lowest pressure in the st 
for P,- Because an objective definition of P, is easily constructed, many 
writers use this value to estimate the intensity of the storm. 

The first approximation of the wind speed within any cyclone is given B 
the gradient wind equation. This can ke expressed in terms similar to Eq 
by: 


Nd pt ae A 5 


g/r + tv, = p(B, - P,) F'(r/R). 


Vg = gradient wind speed 
f =coriolis parameter 
P =density of air 

' = derivative of F with respect to (r/R) 
Near the center of a hurricane Ve2/r is considerably larger than fV,; 
the second term on the left can be neglected within the accuracy of this a 
proximation. Thus the gradient wind, and presumably, the true surface V¥ 
also, is approximately proportional to the square root of (P, - Po). Thi 


lationship is often used in estimating the maximum wind speed in a 
hurricane, (11 
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A hurricane model developed along the lines indicated by Eqs. (1) and (2) is 
described by the following equations: 


P > Py 
apace a agree (3) 
vile + tv, = aaa ire: est (4) 


The symbol R in this model has been identified with the radius of maximum 
wind speed. The agreement between this model and the observations of an ex- 
ceptionally well observed hurricane is shown in Fig. 1. The insert map gives 
the storm track. The dots on the left give the observed pressure at several 
stations in the vicinity of Lake Okeechobee, Florida. The solid line gives the 
theoretical pressure profile fitted to three points within 50 miles of the storm 
center. The corresponding theoretical wind profile is given by the upper curve 
on the right hand side of the figure. The observed winds at one station are 
given by the dots below this curve. A solid line has been drawn through these 
dots by eye in order to obtain a smooth profile. The observed wind speed 
varies in a systematic way from about 65% of the computed wind speed at the 
outer edge of the data shown to almost 90% of the predicted value near the 
zone of maximum wind speed. Equally good agreement between the theoreti- 
cal and observed wind speeds has been obtained in only a few storms. This 
lack of agreement between the theoretical and observed winds is due in part 
to the elementary nature of the model, but perhaps equally to the lack of first 
class wind records near the center of hurricanes. The parameters for this 
model have been computed for a large number of storms and were first 
published in Hydrometeorological Report No. 32.(7) 

Data for new storms are added and data for old storms are corrected as 
new information becomes available. An up-to-date (August 1, 1957) version of 
these data is included in this paper in Table I. 

Eqs. (2) and (4) can be valid only for stationary storms. The wind speeds 
on the right hand side of a moving storm are generally higher than those ob- 
served to the left of the storm track. A simple but incomplete explanation, 
widely accepted for more than half a century, is based on the assumption that 
the forward speed of the storm must be added to the rotary speed on the right 
hand side and subtracted on the left. This would give an asymmetry in the 
wind velocity on the right and left sides of the storm track due to the forward 
speed of the storm. This assumption usually leads to corrections in the right 
direction but not of the proper magnitude. Winds of hurricane intensity are 
significantly affected by friction with the underlying water or land surface. 
The effects of friction and of the inflow necessary to feed the vertical currents, 
prominent throughout the region of hurricane winds, may exceed the effects of 
the forward motion of the storm, at least near the surface. 

Although Eq. (4) is too simple to permit an adequate determination of the 
wind field in a hurricane, it is possible to make use of empirically derived 
relations between the wind speed at nearby locations over water and over land, 
and between the observed wind speed over various types of frictional surfaces 
and the theoretical wind speed to determine a usefully accurate picture of the 
wind field. Figs. 2a and b show the reconstructed pressure and wind fields 
for the New England Hurricane of September 21, 1938.(8) In the derivation of 
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Aug. 28, 1893 28.28 29.61 1.33 958 1003 
Oct. 1, 1893 28.22 29.99 Leer 956 1016 
Oct. 13, 1893 28.33 29.61 1.28 959 1003 
Oct. 2, 1898 28.82 30.03 12k 976 1017 
Oct. 31, 1899 28.70 30.49 1.79 972 1033 
Sept. 8, 1900 27.64 29.78 2.14 936 1009 
Aug. 14, 1901 28.72 30.16 1.44 973 1021 
Sept. 11,1903 28.84 30.12 1.28 977 1020 
June 17, 1906 28.91 29.98 1.07 979 1015 
Sept.17, 1906 28.98 30.38 1.40 981 10293 
Sept.27, 1906 28.50 30.07 1.57 965 1018 
Oct. 18, 1906 28.84 29.80 96 977 1009 
July 21, 1909 28.0) S0ear 1.96 959 1025 
Sept.20, 1909 28.94 30.30 1.36 980 1026 
Cet. 11, 1909 28.30 30.07 het? 958 1018 
Oet.. 17, 1910 27.80 29.19 p RR | 941 989 
"Octe. 18, 2910 28.53 29.77 1.44 959 1008 
Aug. 28, 1911 28.92 30.10 VeL8 979 1019 
Sept. 3, 1913 28.81 29.98 1.27 976 1015 
Aug. 17, 1915 28.01 29.65 1.64 949 1004 
Sept.29, 1915 PARI Te RAO eR aS 2.44 938 1021 
July 5, 1916 28.38 30.03 165 961 1017 
Aug. 18, 1916 cosO0 ose ne Fos a 948 1042 
Oct. 18, 1916 28.76 30.20 1.44 974 1023 
Sept.28, 1917 28.48 29.88 1.40 964 1012 
Sept. 9, 1919 27.44 29.78 2.29 929 1067 
*Sept.14, 1919 Aa near 28.0 inches 
Sept.21, 1920 28.93 29.90 97 980 1013 
June 22, 1921 28.17 30.03 1.86 954 1017 
Oct. 25, 1921 28.29 29.59 1.30 958 1002 
Aug. 25, 1924 28.70 30.33 1.63 972 1027 
*Aug. 26, 1924 28.70 29.62 92 972 1003 
Oct. 19, 1924 28.70 29.82 1.12 972 1010 
*Oct. 20, 1924 28.83 29.62 79 976 1003 
Dec. 2, 1925 28.95 29.90 95 980 1013 
July 28, 1926 28.34 29.91 1.57 960 1013 
Aug. 25, 1926 28.31 30.35 2.04 959 1028 
Sept.18, 1926 27.59 29.99 2.40 934 1016 
*Sept.20, 1926 28620" S0613 1.93 955 1020 
Oct. 20, 1926 27.52 29.97 2.45 932 1015 
Sept.16, 1928 27.62 30.38 2.75 935 1029 
June 28, 1929 28.62 29.97 hess 969 1015 
Sept.28, 1929 28.15 30.08 1.93 953 1019 


Table I - Characteristics of United States Hurricanes 
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STATION 


15 4 28.28 958 Savannah, Ga. 2 
(a 4 28.65 970 Moss Point, Miss. 13 
21 4 28.33 959 South Island, S. C. 1 
zl 4 29.12 986 Jacksonville, Fla. 23 
18 4 28.90 979 Charleston, S.-C. 50 
10 28.48 964 Galveston, Tex. 17 
14 4 29.42(1)996 New Orleans, La. 45 
7 3 29.47 998 Tampa, Fla. 14 
12 4 29.46 998 Jupiter, Fla. 2g 
16 4 29.50 999 Columbia, S.C. 28 
16 4 28.50 965 Ship at Scranton, Miss. 5 
6 4 29.26 991 Jupiter, Fla. 33 
L2 4 29.00 982 Bay City, Tex. 16 
a1 4 29.23 990 New Orleans, La. 43 
10 4 23.36 960 Sand Key, Fla. 7 
11 4 27.80 941 S§.S.Jean Nr. Tortugas, Fla. 0 
ll 4 28.94 980 Tampa, Fla. 45 
- 8 4 29.02 983 Savannah, Ga. 6 
16 4 239.36 994 Raleigh, N.C. 6 
ll 4 28.14 953 Velasco, Tex. 11 
10 4 28.01 949 New Orleans, La.(Pauline St.Wharf) 12 
25 3 28.38 961 Ft. Morgen, Ala. 32 
11 4 28.00 948 Santa Gertrudis, Tex. 6 
2l 4 28.76 974 Pensacola, Fla. ie) 
13 3 28.51 966 Pensacola, Fla. : 12 
8 2 27.44 929 Mean of 2 ships and Dry 
Tortugas, Fla. i 0 
20 4 28.65 970 Corpus Christi, Tex. - 
28 3 28.99 982 Houma, La. 10 
=e! 4 29.37 995 Houston, Tex. 33 
10 4 28.29 958 Tarpon Springs, Fla. 1 
22 4 28.80 975 Hatteras, N. C. 28 
29 4 28.71 972 Nantucket, Mass. 12 
8 4 - - Nr. Dry Tortugas, Fla. (32) (b) 
6 4 29.10 985 Miami, Fla. (0) 
14 4 29.17 988 Wilmington, N. C. 35 
8 4 28.80 975 Merritt Island, Fla. ut 
Loy . 4 28.31 959 Houma, La. 3 
17 4 ZteOk 935 Miami, Fla. 4 
7 4 28.20 955 Perdido Beach, Ala. 1 
16 4 29.16 988 Key West, Fla. 60 
13 4 27.62 935 West Palm Beach, Fla. 3 
15 2 29.12 986 Pt. O'Conner, Tex. 13 
10 4 28.18 954 Long Key, Fla. 7 
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q 
mi 
*Sept.30, 1929 28.80 29.96 1.16 975 1015 40 58 
Aug. 13, 1932 27.83 30.11: 2.28 942 1020 78 12 @ 
Aug. 4, 1933 28.80 29.96 1.16 975 1015 40 24 
Aug. 23, 1933 28.63 29.48 85 970 998 28 54 
Sept. 4, 1933 27.98 29.98 2.00 948 1015 67 13 @ 
Sept. 5, 1933 28.02 30.24 2.22 949 1024 75 30 
Sept.16, 1933 28.25 29.82 1.57 957 1010 53 42 
June 16, 1934 28.52 29.94 1.42 966 1014 48 37 
Sept. 2, 1935 26.35 29.92 3.57 892 1013. 121 6 @ 
"Sept. 4, 1935 28.71 29.89 1.18 972 1012 40 5. 
Nov. 4,1935 28.73 «= - 973 = - - 
July 31, 1936 28.46 30.00 1.54 964 1016 52 19 | 
Sept.18, 1936 28.53 29.42 89 966 996 30 34 
Sept.21, 1938 27.86 29.52 1.66 943 1000 57 50 @ 
Aug. 7, 1940 28.76 29.75 -99 974 1008 34 ll 
Aug. 11, 1940 28.78 30.02 1.24 975 1017 42 26 
Sept.23, 1941 28.31 29.66 1.35 959 1004 45 21 
Oct. 7, 1941 28.98 30.19 1.21 981 1022 41 18 
Aug. 30, 1942 28.07 29.64 1.57 951 1004 53 is @ 
July 27, 1943 28.78 30.02 1.24 975 1017 42 17 @ 
Sept.14, 1944 27.88 30.66 2.78 944 1038 94 49 
*Sept.14, 1944 28.31 29.39 1.08 959 995 36 
Oct. 18, 1944 28.02 29.80 1.78 949 1009 60 
"Oct. 19, 1944 28.42 29.67 1.25 962 1004 42 
Aug. 27, 1945 28.57 30.13 1.56 968 1020 52 
Sept.15, 1945 28.09 30.00 1.91 951 1016 65 
Sept.17, 1947 27.76 292.83 2.08 940 1010 70 
*Sept.19, 1947 28.53 29.70 227 966 1006 40 
Oct. 15, 1947 28.59 29.65 1.06 968 1004 36 
Sept.21, 1948 27.62 29.61 1.99 935 1003 68 
*Sept.22, 1948 28.41 29.83 1.42 962 1010 48 
Oct. 5, 1948 28.85 29.77 92 977. 1008 31 
Aug. 24, 1949 28.86 30.20 1.34 977 1023 46 
Aug. 26, 1949 28.16 30.12 1.96 954 1020 66 
Oct. 3, 1949 28.45 29.95 1.50 963 1014 51 
Aug. 30, 1950 28.92 29.71 79 979 1006 27 
Aug. 30, 1954 28.35 = - 960 + - 
"Auge? Sfp eeG4.  w20.38 fre - 961 - - 
Sept.ll, 1954 27.97 29.77 1.80 947. 1008 61 
Cct. 15, 1954 27.66 29.55 1.89 937 1001 64 
Augs/12) 2656 <28.40. 29277... 2e87 962 1008 46 
Sept.19, 1955 28.51 29.87 1.36 966 1012 46 
Sept.23, 1956 28.76 29.83 1.07 974 1010 36 
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947 
938 
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a ra 
STATION & 
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Panama City, Fla. 

East Columbia, Tex. 
Brownsville, Tex. 

Cape Henry, Va. 
Jupiter, Fla. 
Brownsville, Tex. 
Hatteras, N. C. 
Jeanerette, La. 

Long Key, Fla. 
Columbia, S. C. 

Miami, Fla. 

Ft.Walton, Fla. 

Mean of 2 ships off Hatteras, N.C. 
Hartford, Conn. 

Port Arthur, Tex. 
Savannah WBO, Ga. 
Houston WBO, Tex. 
Carrabelle, Fla. 
Seadrift, Tex. 
Ellington Fld., Tex. 
Hatteras, N.C. au 
Ph. Judith, Re Ls 

Dry Tortugas, Fla. 
Sarasota, Fla. 

Palacios, Tex. 

Homestead, Fla. (FECRR) 
Hillsboro, Fla. 

New Orleans, WBO, La. . 
Savannah, Ga., WBAS 

Key West WBO, Fla. 
Clewiston, Fla. 

Miami, Fla. 

Diamond Shoals L/S, N.C. 
W. Palm Beach WBAS, Fla. 
5 mi. west Freeport, Tex. 
Ft. Morgan, Ala. 

Aircraft Recon. 

Suffolk County AFB, RI- 
40°N = 71°W Aircraft Recon. 
Little River, N.C. 

Ft. Macon, N. C. 

Cherry Point, N. C. } 
Destin, Fla. 
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Note: All values are estimated except pg. 


Po 
Pn 
Pr - Po 
R 


Station - at which pg was observed 
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Footnotes for Table I 


central pressure 

asymptotic pressure 

pressure difference between “outside” of storm and center 
radius to region of maximum wind speed 


(a) computed from pressure profile 
(b) observed from wind speed record 


forward speed of the storm 
time interval over which storm movement is averaged to obtain 


lowest pressure detected by a barometer 


minimum distance from station to track of storm center 


same hurricane as previous line, entering or passing coast at a 
second point 


Lowest in region where pressure profile parameter were com 
29.31 inches observed at Mobile, Ala. later 


Parameters obtained by interpolating is time between ship off 
western end of Cuba and Miami, Fla. and apply to vicinity of D 
Tortugas. 
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Fig. 2. Reconstructed surface pressure, wind speed and deflection ar 
for 1200 EST in the hurricane of September 21, 1938, 
After Myers and Jordan. 
(A) Surface Pressure, (inches of Mercury) 
(B) Wind Speed (M.P.H.). X is location of pressure center 
(C) Deflection Angle (degrees) 
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(4) it is assumed that the winds will blow parallel to the isobars. A more 
histicated theory, as well as observations show that the wind must cross 
isobars toward lower pressure. The angle of inflow, that is the angle with 
ch the wind crosses the isobars, must be asymmetric with respect to the 
ter of a moving storm. The best determination of the angle of inflow for 

7 England hurricanes is shown as Fig. 2c. 8) 


Tides and Storm Surges 


f the effects of wind waves and swell are eliminated, the resulting ele- 
on of the sea surface is usually controlled by the gravitational attraction 
he sun and moon and the topography of the region of observation. Minor 
nges of a seasonal nature are also produced by the seasonal changes in 
ospheric pressure, prevailing winds, and temperature of the water. The 
cts of these gravitational and climatic factors can be computed many 

ths in advance and in this paper are identified as the normal tide. 

n the neighborhood of an intense storm, the actual level of the sea will 

er significantly from the normal tide, and the actual level of the sea may 
-eferred to as a storm tide. The difference between the storm tide and the 
mal tide is called the storm surge. 

fconomically important storm surges may be produced by extratropical 
"ms as well as hurricanes. Extratropical storms have been responsible 
the highest storm tides north of Cape Cod and at a few more southerly lo- 
ons. However, the surge produced by many hurricanes is much higher 

1 those produced by any extratropical storm. Fig. 3 shows the maximum 
orded tide at a representative number of places along the Atlantic and Gulf 
sts of the United States.(9) The remainder of this paper will be devoted 
ropical cyclones and the storm tides produced by them. 

[The peak water level reached in a hurricane will always be somewhat 

ater than the storm tide level because of the effect of waves with periods 
nly a few seconds. It is desirable to separate these two phenomena be- 
se their effects are quite different. If the peak tide level is two feet above 
in sea level, it will penetrate inland to the contour which is approximately 
feet above mean sea level, even if this distance is several miles, and it 
not rise appreciably above two feet regardless of the slope of the beach. 
, two foot wave, on the other hand, may run up to a sloping beach to a 

ht of seven or eight feet above the still water line, but under most con- 
ns the effects of such a wave cannot penetrate inland more than a wave 
th from the still water line, and most of the damage from such waves will 
onfined to a much narrower trip. The waves which accompany the land- 
of a hurricane will generally be in the range of 100 to 1000 feet in length. 
remainder of this paper will be restricted to a discussion of the storm 


n most cases the principal cause of the storm surge is the frictional drag 

yeen the wind and water. Changes in atmospheric pressure, although usu- 
of secondary importance, may occasionally be more important than the 

i effect.(10) Waves breaking against the shore or over barrier reefs may 
contribute to the piling up of water near the shore. 

‘ig. 4 shows the storm surge produced by Hurricane CAROL at 19 Coast 

Geodetic Survey tide stations. The dots on the curves for New London, 

port and Woods Hole indicate approximate surge heights after the tide 
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Fig. 3. High tide data for stations having a record for 20-years or long 


gages became inoperative due to the storm. The dashed line crossing the 
storm surge curves indicates the time at which the storm was nearest e: 
tide station. The storm track is shown to the left of the storm surge cur 
The letters on the map and beside the curves identify the storm surge cul 
with the geographical location of the tide gage. The maximum surge he 
and mean tide range are shown on the left. 

With the exception of Willets Point, near the western end of Long Is 
Sound, the maximum surge height at all stations between Hampton Roads 
Cape Cod occurred between the passage of the storm and two hours later, 
Willets Point, the maximum surge occurred more than four hours after tf 
passage of the storm center. Thus, the water continued to rise at many 
to the left of the storm track after the wind had acquired an offshore cor 
ponent, suggesting that the surge arrives as a wave which had been genet 
elsewhere. This is particularly true at Willets Point. The notion of a W 
receives further support from the damped oscillations readily seen in th 
records for Atlantic City, Sandy Hook, and the Battery, and suggested in 
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ig. 4. Track of hurricane and storm surge curves for hurricane CAROL, 
1954, 


yrds of several other stations. These oscillations, called resurgences by 
field,(4) are characteristic of the hurricane surge in this region. 

‘ig. 5 shows the storm surge as a function of time at the recording tide 
ion nearest the point of maximum storm surge height for several hurri- 
2s. A trend toward rising or falling sea level may become evident 12 to 
iours before the arrival of the storm. In either, a rapid rise begins when 
storm is about 100 to 200 miles distant, that is when the wind speed reach- 
bout 30 miles per hour with the approach of the storm. This period of 

d rise lasts from two to six hours and is generally followed by a slightly 
e rapid drop to something below the normal water level in areas with good 
nage. In marshland or other areas with poor drainage, the water may re- 
n well above normal for several days. Fig. 6 shows the contours of the 

m surge along the coast for four hurricanes. Notice the tendency for the 

t value to be higher on the right or to extend farther from the center on 
right hand side of the storm track. Notice also that the peak storm tide 
‘be as much as 10 feet higher than the storm tide 60 miles to either side 
le peak, 
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Fig. 7 based on data furnished by the Corps of Engineers, Coast and G 
detic Survey, Woods Hole Oceanographic Institution and others, shows the ~ 
storm surge heights produced by Hurricane CAROL at a great many locatior 
along the southern New England and Long Island coasts. Here the effects 
topography and local wind conditions are shown by the higher surge heights 
near the heads of converging bays and at exposed points, and the decreas 
surge heights near the widest section of Long Island Sound. The varied 
of topography and local wind conditions make it difficult to construct eraplis 
such as Fig. 6, or to develop any reiation between storm surge height and 
meteorological parameters. It is clear that one should avoid using surge ; 
heights obtained near the heads of estuaries in developing such relationship: 
but it is not possible to avoid all local effects. & 


Maximum Storm Tide as a Function of Storm Intensity | 


The effect of wind in piling up water is nearly proportional to the wind } 
stress. The wind stress is generally taken as being proportional to the 
of the wind speed, and according to Eqs. (2) and (4), this is approximately 
portional to the pressure deficiency in a hurricane. The direct hydrosta 
fect of the decreased pressure is also proportional to the decrease. Thus, ? 
are encouraged to look for a correlation between (P, - Po), or Po alone, 
the maximum storm surge, that is the maximum effect of the storm on re 
level of the sea. The results of such an effort are shown in Fig. 8. 4 

The highest observed tides along the Gulf Coast of the United States duri 
30 hurricanes is indicated by dots as a function of the central pressure of t 
storm as it crossed the coast. An attempt has not been made to eliminatet 
normal tide, as the normal tide range is less than two feet in most of thigl 
region, and sufficient data for the elimination of the normal tide were not i 
available for all of the reports. The regression equation based on these 
gives a correlation coefficient of 0.68. The extreme differences between 
storm and the normal tide for 10 Atlantic Coast hurricanes are indicated bj 
x’s. The tide range along the Atlantic Coast is much greater than in the Gul 
of Mexico, and a correction for the normal tide is necessary in order to ok 
tain any useful degree of correlation. Largely because the variability duet 
the normal tide is eliminated, the correlation coefficient for Atlantic data! 
been increased to 0.72 even though these data were not considered in the dé 
vation. Consideration of this single parameter is thus sufficient to account 
for half of the observed variability in maximum storm surge heights, (12) I 
improvement was obtained by considering P, and R as defined in Eqs. (3) a 
(4), or by using the storm tide potential as defined by Reid. (6) 


? 


Effects of the Forward Motion of the Storm 
Proudman(13) using a very elementary theory, has shown that the respe 


of the water in a rectangular canal to a moving atmospheric disturbance ¢ 
be expressed in the form: 
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h_ = height of water level disturbance 


h = Atmospheric pressure disturbance expressed as the height of a 
water barometer 


V, = speed of atmospheric disturbance 
D_ = depth of canal 
g =acceleration of gravity 


ter writers(14) have shown that these results are modified somewhat, when 
e effects of the earth’s rotation and the two dimensional nature of the dis- 
rbance are considered. The essential feature of the development which 
rries over to more complex situations is, that for any basin there is a criti- 
J velocity of motion of the atmospheric disturbance which will produce maxi- 
im response in the water. An examination of the residuals obtained from 

g. 8 supports this concept, but there are not sufficient data from any type of 
astline to establish the proper form of this relationship from empirical 
idence. 

The speed of the storm can become important in quite another way, for the 
eed of the storm, together with the radius of maximum winds determines 

= duration time and the length of the fetch over which the storm winds can 

t to produce storm tides. 


Hurricane Frequency 


From a practical point of view, it is necessary to consider the frequency of 
rms as a function of intensity, size, and region of occurrence. Any real 
ends in such frequencies would also be highly significant in an economic 
nse. Unfortunately, the amount of data available for making such studies is 
ther limited and all data are not of the same quality. The nature and sense 
the available data can be indicated by three figures adapted from recent 
sather Bureau reports to Congress. (16) 

Fig. 9 shows the frequencies of tropical storms and hurricanes for the 
riod 1887-1956. A tendency for years of above normal hurricane frequency 
be followed by a period of below normal hurricane frequency can be easily 
ted. There is no evidence of any well marked trend or regular cycle in 
rricane frequency. It may be possible to obtain a useful indication of the 
mber of storms to be expected during the next two or three years by con- 
lering the number which have occurred during the last 10 to 20 years, but 
» best estimate of hurricane probabilities more than a very few years in the 
ure would be obtained from considering the frequencies for the entire peri- 
of reliable records. 

Fig. 10 shows the number of times destruction has been caused by tropical 
rms on the United States mainland during the present century. 

Fig. 11 shows the accumulated frequencies of tropical storms, whose 
itral pressures were 29.00 inches of mercury or less, for several sections 
the U. S. Coast. Cape Hatteras is taken as the dividing line between the 
rth Atlantic and South Atlantic sections of the coast for this figure. The 
a for the Florida Keys are clearly different from those of the other regions. 
ither the size of the sample available, nor the present state of theoretical 
ywledge concerning hurricanes is sufficient to determine whether the 
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Chronology of Tropical Storm Tracks 
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Fig. 9. Chronology of storm tracks. The occurrence of hurricanes a 
tropical storms in the different coastal regions shown at 
left, and the total occurrence. 
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Fig. 10. Number of times destruction was caused by tropical storms, 
1901-1955. 


fferences in the other five curves are due to sampling error or to some 
ndamental aspect of hurricane behavior. 


Hurricane Observations and Forecasts 


If a hurricane approaches land in the neighborhood of a powerful weather 
dar installation, such as that now operating at Cape Hatteras, N. C. and on 
der for about 15 other coastal sites, an experienced radar hurricane ob- 
rver can locate the center of the storm within 5 to 10 miles while the storm 
nter is-within 150 miles of the station. The accuracy decreases as the 
stance from the radar station is increased. 

When the storm is out of range of land based radar, the position of the 
rm center must be determined from the available ship and aircraft reports. 
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casionally the eye of the storm will pass directly over a ship. Somewhat 
re frequently, but not on every day of the hurricane’s existence, an airplane 
1 penetrate the eye of the hurricane. In these cases the storm center can 
en be located within an accuracy of about 10 miles. In the more usual case, 
storm center must be located in the center of an isobaric system de- 
‘mined from fewer than a dozen ship reports, or from radar observations 
de from aircraft on the edge of the storm. The error in locating the storm 
iter in this way will frequently be as much as 20 or 30 miles and occasion- 
y as much as 60 miles. Errors in communication via radio and teletype- 
iter compound the uncertainties resulting from the currently unavoidable 
yrtcomings in our observational procedures. 
A few months after the storm, when mailed in reports, free of communi- 
ion errors, as well as much additional information not available to the fore- 
ster during the storm are obtained, it is possible to repeat the analysis and 
ermine the position of the storm track with an accuracy impossible during 
+ storm. 
A simple extrapolation of the storm path accounts for considerably more 
n half of the obtainable skill in hurricane forecasting. Even the more 
Ihisticated forecast methods, based on a knowledge of flow patterns at all 
els, give only a velocity and acceleration of the storm movement, and per- 
uates any existing error in the storm position. 
When all of these factors are considered, it is evident that storm 
vements cannot be forecast with timetable accuracy. A study of forecast 
suracy for the period 1952 to 1955 indicated an average error of 69 nautical 
les in the 12-hour forecast of storm position, (1 5) This was made up of an 
srage error of 83 miles during the period 1952-1953 and 63 miles in 1954- 
95. A gradual improvement in forecast accuracy is to be expected as the 
sult of improved instrument and forecast techniques. 
A hurricane watch is announced for specific areas whenever it is recognized 
t a specific hurricane threat exists for any region of the coast. The hurri- 
1e watch is announced before it is possible to delineate the area in which 
-storm is likely to go inland and when the hurricane effects are believed to 
at least 24 hours away from the coast. A hurricane warning is issued for 
pecific area whenever there exists a reasonable probability that the storm 
1 go inland in that region within the next 24 hours, When the average accu- 
-y of forecasts and the length of coastline seriously affected by a hurricane 
> considered, it is evident that the hurricane warnings must cover an area 
eral times as large as that which will actually suffer severe damage from 
storm. 


SUMMARY 


A full realization of the hurricane’s potential for danger and of the difficulty 
redicting its future behavior exactly, should materially reduce the loss of 
»and property due to these storms. 

More than 390(1) persons lost their lives in Hurricane AUDREY, largely 
ause they had lived through other hurricanes with comparatively little 

nage and did not believe that this storm could be as severe as predicted. 
Hurricanes vary in intensity. Tropical storms which barely qualify as 
ricanes may elevate the sea level no more than four feet, but severe hurri- 
es may elevate the sea level by 15 feet or more above normal tide for the 
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time of the hurricane landfall. Intense extratropical storms may product 
rises of five to six feet above normal. 

The sea level may be raised by one or two feet for a thousand miles of 
coastline by even a moderately intense hurricane, but the really severe sto 
surges of more than 5 feet above the normal tide level rarely extend oval 
more than 100 miles of coastline. 

‘The duration of the rapid rise in sea level is rarely greater than 6 hour: 
and in areas with good drainage the fall is even more rapid. In areas of 
drainage, moderately high water may persist for days after the storm. : 

Knowledge of the intensity of a hurricane before it reaches the coast 
meager. The average error in a 12-hour forecast is approximately 60 naut 
cal miles, and if warnings are to be effective, an area several times as lar 
as that actually affected must be alerted. 

Disasters such as that which occurred with Hurricane AUDREY can 
eliminated by the evacuation of the coastal plain in the threatened region & 
elevation of 10 to 20 feet above the normal high water level. In marshla 
where high ground may be as much as 40 to 100 miles away by highway, m2 
mass evacuations may be required to avoid one catastrophe. With inp 
instrumentation and forecasting techniques, the number of evacuations, 
proved to be unnecessary, will be reduced, but no technique now avaite 
permit the prediction of hurricane motion with timetable accuracy. 
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MECHANICS OF SEDIMENT-RIPPLE FORMATION®2 


Corrections by J. Bogardi 


. BOGARDI.!—Corrections to Discussion by J. Bogardi: 


age 1832-3, 14th line “. . . difficulties to Mr. Liu—of the multitude...” 
iid be “difficulties to Mr. Liu to determine these correlations—of the 
titude.” 


lage 1832-4, Eq. (8) “T, = U*2 = ” should be “T, = pU*= ”. 


forrections to Closure by H. K. Liu: 


age 1832-29, 16th line “. . . This factor is proportional to gd/U, ane ee 
iid be “The third parainciey B is proportional to d/DS = gd/U, 2, hie! 
inel stability factor, which is the inverse...”. 


age 1832-29, delete (a) the words “in the form” before Eq.(35), (b) E 

and (c) the line following Eq. (35), and substitute the following: 

If one uses the channel stability factor d/DS = gd/U, 2 the configuration of 
ged is quantitatively defined by ne parameter B mentioned above unique- 
The relationship of 8 and gd/U,?2 according to Bogardi is 


pe 2a 
= a 

Y Wd 
rom the data of Fig. A, Bogardi assumes that N = 0,882. 
ogardi further states: Obviously there is a parameter o@ proportional to 
V , which defines the bed configuration uniquely too and 


Ss 
4 


(35) 


(36) 


inally Bogardi pointed out that a third parameter € proportional to U, 
rally defines the bed configuration uniquely too: 

/ 

Pi (37) 

fea 
os 
, @ and € are obviously interrelated, and their value depends upon the 
erature, as well as on the specific gravity of the sediment. Assuming a 
erature of 20° C., and a specific gravity of 2,65, the parameter values 
ybtainable directly from Eqs. (35), (36), and (37). 

1 


ig. 19 shows after Bogardi the U, = dey, equations. 


roc. Paper 1197, April, 1957, by H. K. Liu. 
sst. Prof. Technical Univ. of Budapest, Hungary. 
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According to Fig. 19 ripples are only possible if 1,777< € < 4,06 an : 
d< 2mm. Bogardi constructed similar figures for Eqs. (35) and (36). 
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NORTHEASTERN FLOODS OF 1955: RAINFALL AND RUNOFF2 


Discussion by Howard M. Turner 


TOWARD M. TURNER, ! M. ASCE.—This is a very interesting engineering 
cription of the floods of August and October, 1955. The data selected and 
presentation of it give an excellent summary of the magnitude and extent 
his storm, which, as the author points out, exceeded previous records of 
nsity of rainfall and flood flows and covered a tremendously wide area. 

- writer has made some further comparisons of the magnitude of this flood 
1 other large floods in New England. 

Fig. 1 shows area-depth curves for total rainfall in various New England 
"ms showing how the 1955 storm exceeded previous records even of the 
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very large four-day hurricane storm from September 17-21, 1938. Most of 
the 1955 storm came the two days of August 18 and 19. j 
The Storm Rainfall report of the Miami Conservancy District, Revise d, 


the country from 1869 to 1933. On Fig. 2 the 1938 and 1955 hurricane sto 7 
are plotted on the two-day area-depth curves for the Northeast taken fron 4 
that report showing how much the 1955 storm exceeded any of these recor¢ 
storms. 


Fig. 3 shows a curve of high rainfall frequency taken from “The Frequency 
of High Rates of Rainfall” by Allen Hazen, published in 1921.(2) The writ 
has plotted on this David L. Yarnell’s values, published in 1935(3) for 10-_ 
year and 100-year maxima estimated for Massachusetts and also the maxi- 
mum actual rainfalls for different periods of time recorded during the 195 
storm. The results of the latter would require an entire recasting of these 
frequencies derived thirty-six years ago. 
In a paper of “Flood Flow Formulas based on Drainage Basin Character= 
istics” by H. B. Kinnison and B. R. Colby, (4) a frequency curve of storm r2 
fall was included. This is shown on Fig. 4. On this are plotted the maximi 
records for the 1955 storm showing the extreme magnitude of this stormé 
pared with that curve. : , 
It is naturally expected that a longer period of records would have the r 
sult of increasing the maximum results obtained, but it is interesting to no 
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the present period does not seem to show any orderly increase as a long- 
ime of records alone would cause. Fig. 5 shows records of large three- 
storm rainfall since 1855. The two largest storms have come in the last 
ears, exceeding by 25 per cent anything of record before. 

yonsidering flood flows, records of the gage heights on the lower 

necticut River at Hartford are available for the last 114 years. There are 
a few high-gage records going back for 318 years. The greatest flood 

in 1936; the next, in 1938; and now the third, in 1955. As far as can be 
rmined from storm and flood histories, nothing of these recent large 
nitudes occurred in over 300 years. Fig. 6 shows this large flood record 
ie Connecticut River at Hartford back to 1683. From all that can be 
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srmined, the three largest floods have occurred in the last twenty years. 
great size of the 1955 storm is shown by the large flood it caused on this 
sr though it was largely confined to the lower third of the drainage area. 
[fhe author has shown the relation between the peak flow and drainage 

as on his Table 5 and Fig. 6. The writer has made up Fig. 7 showing the 
same data for New England, including the large floods previous to that of 
ust, 1955. A curve of the 40% Myers rating is shown as an enveloping 

ye including practically all records up to that time. The 1955 flood re- 
es a 70% rating curve to include the high records of the 1955 flood for the 
nage areas of 300 to 700 square miles. 
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These very high rates of flow on these fairly small drainage areas sho 
in the Myers rating and in the rainfall frequency studies are due to the he y 
rainfall in a short time. The author shows on Fig. 3 the rainfall curve at 


ally 97 per cent of it fell in 28 hours; and about 83 per cent of it, in 24 hou r 


Runoff 


The author gives some figures of the percentage of runoff to rainfall, 
which, except for a few cases, appear to be somewhat lower than might be 
pected from a flood of this intensity. Some of the very high ones that are © 
given may very likely have been affected by the failure of dams. 

The author has noted the rapid rise and fall of the hydrographs of this 
flood. In the case of the large rivers the reason for this is quite evident. 
For example, on the Connecticut River the heavy storm flows were con C 
to the lower third of the drainage area. Channel storage was thus confined 
only the lower portion of the river valley, and the recession curve of the 
hydrograph representing the drainage of the channel storage would be short 
in duration as the amount of storage to drain away would not be so great. 1 
same thing shows in the recession limb of the hydrograph of the Delaware 
River at Port Jervis, where the heavy storm rainfall was concentrated into 
about 50 per cent of the drainage area. 

The hydrographs of the smaller rivers remain to be analyzed. The wri 
believes, however, that the 6- or 12-hour unit hydrographs for these rivers 
will prove to be close to the corresponding unit hydrographs for the 1938 
flood except that the peaks may be somewhat higher. 


Frequency 


The question of frequency is one of the most difficult problems in flooc 
studies, but some frequency rating of the size of a given flood seems to be 
quired. The author gives a tabulation and a frequency curve for the No , 
east, according to the United States Geological Survey’s method of rating t 
recurrence interval of any given flood as its ratio to the mean annual floot 
but the curve gives the recurrence interval up to 100 years only, not long 
enough for many requirements. This curve shows very different results f 
that developed for Connecticut by B. L. Bigwood and M. P. Thomas. (5) 

The section on flood frequency of the 1942 report of the Committee on 
Floods of the Boston Society of Civil Engineers, published in 1942, (6) analy 
the frequency of floods on the Connecticut River at Thompsonville, where 
there were then over 93 years of record. Various methods were used. T 
results, using data up to 1935, i.e., prior to the 1936 and 1938 floods, gave 
frequency of the 1936 flood on 10,000 to 16,7000 years, depending on th 
method used. The same computations made with three years’ more data 
cluding these two large floods gave frequencies from 300 to 400 years. 

Apparently the same difficulty is still found by the United States Corps 
Engineers as shown in its recent report on the 1955 flood, (7) where for fc 
rivers in southern New England the flow for a given frequency used in 19% 
was much higher than shown by the curves tentatively adopted prior to the 
1955 flood. See Table 1. 

The paper by Kinnison and Colby, published in 1945, (4) which has been 
a great deal in Massachusetts, rates floods as “minor, ‘2 “major,” and “ré 
with stated frequencies of 15, 100, and 1,000 years respectively, On the 
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on which there are records concerning the 1955 flood, the Kinnison-Colb; 
figures would be as shown on Table 1. These generally rate the higher 1 
flood flows as equal to or exceeding a “rare” flood. 

In that paper the flood frequency ratings were based upon precipitatior 
depth of runoff. The use of precipitation in determining flood frequency 
its difficulties because of the varying percentage that occurs in the form 
runoff in the flood, but the rainfall data are more inclusive. The writer 
lieves that the flood runoff figures would have some advantage as data fo 
flood frequency computations instead of the peak floods, The runoff data; 
often not so readily available and require more work to use, but the resu 
are on a broader basis, eliminating the varying characteristics of each 
stream. With a runoff frequency established, the peak flood of each locatik 
could be determined by unit hydrograph methods. On a few of the streams 
for which Kinnson-Colby figures are available for the high 1955 flood area 
the runoff compares to their “rare” or 1,000-year flood. 

All this points out the difficulty of the problem... Each of these success 
floods gives more information with which to work; but as the author points 
out in his paper, figures that still show such wide variation are not very 
“satisfactory.” Yet, it is a problem which is important. In over 950 yea 
of records of the flow of the river Danube, the greatest flood was in 1501; 
next in 1787 was 85 per cent of it; and the third one in 1899, 75 per cent 4 
it.) There is no question but what these extraordinary events occur at’ 
intervals, and we have to know something of the degree, time, and magnit 
in flood works. 

Fortunately, frequency is not required so much in the design for spill 
capacities of dams except indirectly; but for other designs, such as flood: 
control reservoirs, flood channel improvements, and other flood-control 
work, some method is needed to fit the design flood to the value of the pr 
posed construction. The 1955 flood has certainly given perhaps an uns¢c! 
tific but in some ways very convenient yardstick. It has been taken as th 
design flood on some recent river-improvement designs. As far as the | 
is concerned, a large storm like that of 1955 is, of course, a very good ff 
measure, It is much easier to describe flood-control measures for a fle 
intensity that has already been experienced than it is to talk about the so 
what mythical frequency of a 100- or 1,000-year flood, which may come: 
next year. 

It is interesting to look back and see the tremendous development of # 
science of flood analysis and prediction since the 1927 flood. There are 
problems, but the progress made justifies the belief that the future will s 
them. 


REFERENCES 


1, “Storm Rainfall of Eastern United States”, State of Ohio—Miami Con- 
servancy District, Technical Reports—Part V—1936. 


2. “Frequency of High Rates of Rainfall”, by Allen Hazen, Engineering q 
Record, Vol. 87, No. 21, 1921. 


3. “Rainfall Intensity- Frequency Data”, by David L. Yarnell, U. S. Dept 
Agriculture—Miscellaneous Publication—No. 204. 


-E DISCUSSION 1880-15 


“Flood Formula Based on Drainage Basin Characteristics”, by H. B. 
Kinnison and B. R. Colby, ASCE Trans. Vol. 110 (1945), p. 849. 


“A Flood Flow Formula for Connecticut”, by B. L. Bigwood and M. P. 
Thomas, U. S. Geological Survey Circular No. 365, Washington, D. C.— 
1955. 


“Journal of the Boston Society of Civil Engineers”, Report of the Com- 
mittee on Floods, Volume XXIX—No. 1—Section 2. 


“New England Floods of 1955”, Part 2—Flood Discharges, Corps of Engi- 
neers, U. S. Army—Office of the Division Engineer, New England Division 
—Boston, Massachusetts, April, 1956. 


2r 1890 December, 1958 HY 7 


ENED 
Journal of the 


HYDRAULICS DIVISION 


Proceedings of the American Society of Civil Engineers 
a 
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SYNOPSIS 


leasurements of distribution of air concentration in self-aerated flows 
presented in this paper. The experiments were made in a rough channel 
ind-grain type surface at various slopes and discharges, and the data 

2 used as a basis for study of the mechanism of entrainment of air and to 
te the air content and distribution to the characteristics of the flow. The 
ysis of the data shows that the air distribution can be adequately de- 

bed by relationships based upon a simplified concept of turbulent transport 
thus are functions of the flow characteristics. The maximum depth and 
nean velocity are both shown to increase above those of a corresponding 
erated flow. 


INTRODUCTION 


characteristic of high-velocity, open-channel flow is the phenomenon of 
aeration in which atmospheric air is insufflated into and mixed with the 
to create the appearance of “white water” with its violently agitated and 
efined free surface. This condition is frequently observed in flows down 
) chutes and spillways. 

has been reasonably well established by several observers(1,2,3) that 
lows over a spillway from a quiescent reservoir, incipient aeration does 
ccur on the slope until a point or region is reached at which the boundary 
> thickness is equal to the depth. The characteristic roughening of the 
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water surface immediately upstream of the appearance of “white water” 
be readily observed in such occurrences, With increasing initial water 
depths, the points where the surface roughening and the “white water” occ 
move downstream along with the intersection of the boundary layer with t 
water surface. These considerations suggest that the aeration phenomenz 
related to the conditions of turbulence in the flow. Examination of such ik 
by means of high-speed photography indicates that the upper boundary of 
flow is rather ill-defined. It consists of a zone which appears white beca 
of entrained air, while above this zone a spray of water droplets occurs mic 
ing more or less parallel to the flow and below which is a region of discre 
air bubbles suspended in the fluid. Other investigators have described air 
entrained flow as being made up of air bubbles entrained in water in the lh 
levels and water particles in the air in the upper regions. This concept 
further emphasized by measurements of the air concentration distri- 
bution(4,5,6) which show that the air concentration increases continuous y 
from the bed, exhibiting a smooth transition from a finite value near the be 
to a maximum value of 100 per cent at the free surface which value is ap 
proached asymptotically. 

This paper presents the results of systematic experiments on distribu 
of air in self-aerated flows in a channel set at various slopes and opera 
various discharges. A provisional analysis is made of the data with resp 
to air distribution in a vertical section transverse to the direction of floy 
with a view of arriving at empirical relationships between concentration 
parameters and the flow characteristics. 


O 


Experiments and Results 


Apparatus 


The experiments were undertaken in a channel 50 ft. long, 1.5 ft wide 
1 ft deep, which could be adjusted to any slope from horizontal to nearly 
cal. The width was chosen so that for ordinary flows sidewall effects due 
growth of the sidewall boundary layer and air entrainment along the sides 
could be avoided in the central portion of the cross section. That the wic 
was adequate to make the flow two-dimensional in the center portion was 
demonstrated by earlier measurements of air concentration and velocity 
transverse section at the end of the flume.(4) The flow was obtained by g 
ty from the laboratory main supply flume through two feed lines and cont 
by hydraulically operated valves. The water reached the inlet of the flun 
through hollow support members and through a rectangular conduit on the 
underside of the channel. Fig. 1 shows the channel and appurtenances. 

At the inlet the flow is guided through two 90-degree vaned bends and 
contraction. The inlet gate itself is cantilevered from its mounting at th 
of the flume so that there are no slots or guides to disturb the flow. Ful 
guidance is provided by a lip attached to the bottom of the gate and exter 
inward to further contract the flow; therefore, the gate opening correspc 
to the depth of flow at the upstream end of the channel. The maximum ¢ 
opening is 0.5 ft. Pitot measurements of the flow just downstream of th 
with an 0.25-ft gate opening showed no deviation greater than 2 per cent 
the mean velocity except in the bottom corners where slightly smaller vi 
ties were encountered. “<4 

In the experiments here described, the air-entrainment process was i 
tensified by the installation of artificial roughness on the channel bed thé 
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siderably rougher than the painted steel sidewalls of the channel itself. 

material used for the roughness was a commercial non-slip fabric coat- 

vith granular particles. The particles, whose mean size was 0.028 in., 

a mean spacing of 0.039 in. and were embedded in a mastic, as shown in 

m2. 

The discharge rate was measured by two calibrated Pitot cylinders mount- 

n the supply lines and leading to differential gages mounted on the control 

el. Each discharge meter was calibrated separately for flows from 1/2 

0 cfs with an accuracy of the order of 1 per cent. 

\ control panel contained the actuating components for the hydraulic 

tem and indicators for the measurement of slope, gate opening, and dis- 

rge which permitted adjustment of these variables to predetermined 

les. 

The electrical air-concentration measuring instrument has been fully de- 

ibed in an earlier report. (7) The method consists basically of a measure- 

it of the difference between the conductivity of a mixture of air and water 
the conductivity of the ambient water alone. A strut supporting a pair of 

bes is combined with the electrical circuit so that the air-concentration 

isurements may be made in a small region of the flow. The probes were 
in. wide and 1/4 in. apart. The instrument is direct reading, and the 

ition between conductivity and air concentration can be determined ana- 

cally. The strut holding the probes is arranged so that it is possible to 

rerse the flow cross section both vertically and laterally. The instrument 

imited to measurements not closer than 0.02 ft from the walls because of - 

rent deflections by the walls at lesser distances. Since the measurements 

ne air concentration depend upon the resistance or specific sensitivity of 

suspending water, the meter was zeroed in water continually circulated 

n the test channel through a deaerating tank, thus providing water of the 

ie character as that in the flume. Oscilloscope measurements show a 

idly fluctuating current between the probes. Since the meter could not 

yw these rapid fluctuations, the value read on the meter represented an 

rage value. This average compared well with results obtained by means 
mechanical sampler which had been previously developed at the Labora- 

‘and which was not subject to the same averaging problems. A photograph 

1e probes and strut supports for the air-concentration meter is shown in 
3. 

ome of the experimental runs also required taking velocity traverses in 

icals through the aerated flows. The instrumentation for making these 

srvations has been previously described(8) in some detail. It consists 

ntially of a device for timing the travel of minute salt-water cloudlets in- 

ed into the air-water mixture. The operating principle of the velocity 

er consists basically of marking a small element of flowing mixture and 

recording the time interval required for this marked element to traverse 

‘ed distance. The marking is accomplished by making a diminutive portion 

owing water more highly conductive to electrical current by injecting a 

li amount of salt solution into the stream. The passage of this ionized 

diet is then detected by electrodes at stations fixed 3 in. apart in the flow 

. The injections were made at a speed of about 20 per sec with individual 

stions of about 0.03 cu cm of 6-per cent saline solution injected in about 

(0 sec. The accuracy of the air-concentration meter and the velocity 

2r was checked by comparing the measured water rate entering the 
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channel with the integrated flow determined by the velocity traverses 2 
air-concentration traverses over the channel cross section. 


Experimental Data 


For the air-distribution measurements described in this report, the 
channel slope and water discharge were varied independently. Experiments 
were made for slopes of 7.5, 15, 22.5, 30, 37.5, 45, 60, and 75 degrees. Fo 
each slope, measurements of air concentration were made for total-water 
discharges from 2.2 to 9.6 cis and for some slopes up to 15.0 cfs. " 

The measurements were made at Station 45 (45 ft from the inlet) in a ve 
cal centerline plane. Readings were taken at 0.01-ft intervals on a line 
normal to the channel bottom with the lowest point 0.02 ft from the bed 2 
the uppermost point where the probes were completely out of the flow andt 
meter registered zero water concentration. The concentration profile was 
obtained for each slope and each water discharge after the channel had be 
adjusted for equilibrium conditions. For this purpose equilibrium flow was 
considered to be that condition of the flow of the air-water mixture for whit 
the air-concentration distribution was the same at two sections 10 ft apar 
along the channel. This condition could be obtained by the adjustment of 
initial depth of flow through control of the inlet-gate opening. Repeated 
measurements of the profiles were made at the two sections for different 
openings until the two air-concentration profiles were similar. Compariso 
were made of the concentration at corresponding distances normal to the 
A typical comparison of the profiles at the two stations for equilibrium f 
is shown in Fig. 4. 

The results of the experiments in the form of measured concentrations 
each slope and discharge have been tabulated in Appendix B for reference. 
tabulation of computed parameters is given in Appendix C. 

The air concentration C is defined as the ratio of the volume of air per 
unit volume of air and water as measured by the concentration-meter prol 
It is assumed that the concentration as measured represents the average - 
value of the concentration in the region described by the probes of the con 
tration meter and is applicable to the midpoint of the probes. Essentially 
continuous curves of air concentration with respect to distance normal to 
bed were obtained for each experiment. The concentration increases gra 
ly from the bed—more rapidly in the central portion and then more slowly 
the upper region, apparently asymptotically approaching the limit of 100 p 
cent air concentration, as is shown in Figs. 4, 8, and 9. It appears that tl 
complete curve is comprised of two parts which have basically different 
characteristics. Although there is no sharp demarcation between the two 
parts, the curves tend to support the description that in the lower regio s 
the stream air bubbles are suspended in water, while the upper regions 
sist of water droplets in air. The existence of these water droplets can1 
ly be detected by holding one’s hand just above the main flow. 

Since there is no definite upper boundary to the air-entrained flow in 
open channel, questions of the definition of the linear dimensions of the f 
arise. A number of different depth parameters related to the air concent: 
and distribution may be defined, each having particular applications in de 
scribing the bulk-flow conditions. These depth parameters may be descr 
as follows: 


CE OPEN CHANNELS 1890-5 


(1) A depth d is used to represent a mean depth of flow that would exist if 
all of the entrained air were removed up to the highest point where 
water is found. It therefore corresponds to the depth of a nonaerated 
flow of a given discharge with a velocity equal to that of the entrained 
flow. It may be defined as 


d = [(I-C)dy (1) 


where y is the normal distance from the bed and C is the concentration 
of air as a function of y. The value of d can be approximated with suf- 
ficient accuracy by a summation of the products of the water concen- 
trations and their corresponding incremental depths. 

(2) An upper boundary of the air-entrained flow or an upper depth d, may 
be defined as the value of y where the air concentration has some arbi- 
trarily prescribed value. For some purposes, such as determining the 
mean depth, and because of the asymptotic nature of the curve, itis | 
convenient to define dy as that value of y where the air concentration C 
is 0.99. Other depths may be similarly defined. Such a depth has 
previously been defined ( ( as that corresponding to that value of y 
where the air concentration is C = 0.95. This depth was chosen arbi- 
trarily as a depth where the concentration of air can be reliably 
measured and which includes from 98 to 99 per cent of all of the water 
and consequently represents a maximum distance from the bed where 
water will be found. 

(3) A third depth parameter may be defined on the basis of the air- 
concentration distribution as that depth dp which represents the value 
of y where the transition from the distribution in the lower regions to 
that in the upper regions occurs. Its value would then depend upon the 
characteristics of the channel as well as the air distribution. This 
depth parameter will be discussed in detail in a later section. 


The mean air concentration in the vertical can also be defined in different 
ys depending upon the particular application. The mean concentration over 
- whole range of concentrations measured is defined as 


= l U 
2 ig comeing BAKE (2) 
dy i y 


ere C is the mean concentration of the entire flow and dy is the so-called 
yer limit where the concentration is 0.99. Another mean concentration Cy 
y be defined as the mean concentration in the region below the transition 
th y = dy which applies to that air which is being transported by the flow, 


tis 
d, 


s | (3) 
C, = all Cdy 
TO 
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Dimensional plots of d, d,, and dy for the various discharges and for ¢ 
slope are presented in Fig. 5 to indicate the range of values as determineg 
from the concentration data. In Fig. 5(c) the mean depth d varies with dis 
charge and slope in much the same manner as nonaerated flow. The meg 
depth increases regularly with the discharge and decreases with the slope. 
The data for the upper limit dy (Fig. 5(a)) indicate that for a given discharge 
d,, tends to decrease at first and then to increase as the slope increases; ~ 
consequently, the maximum depth tends to be nearly the same for correspon 
ing discharge rates on the various lower slopes. It appears that the incre; 
depth due to aeration increases with both slope and discharge; however, 
combination with the mean depth, which decreases with slope and increase 
with discharge, the upper limit tends to decrease and then increase as the 
slope increases. As a result, over a considerable range of slopes (15 to 3 
degrees), the upper limit is essentially the same and depends only on dis- 
charge. A similar situation exists in regard to dp (Fig. 5(b)) in that it als 
shows a tendency to decrease gradually and then to increase gradually as 
slope increases. Consequently, the range of variation of dp for any given 
charge is relatively little over a wide range in slopes. 

The mean concentration C for the entire flow and Cy, the mean concen= 
tration in the region below dy, have been plotted in Figs. 6(c) and 6(b) in 
terms of the discharge and for each slope. It is interesting to note that the 
mean concentration for the entire flow is nearly constant for a quite wide 
range of discharges but increases very considerably with the slope. The 
mean concentration Cy for the region below dy has similar characteristic 
In Fig. 6(a) the concentration Cp at the transition depth is plotted. This 
concentration varies but little with the depth but increases with the channel 
slope. 


Analysis of the Data 


Distribution of Air in Aerated Flow 


A qualitative concept of air-entrainment phenomena has been establisht 
with the aid of the air-concentration determinations made transverse to th 
direction of self-aerated flows, high-speed photographs and motion pictur 
of the flow, and other similar information. Since it has been shown that aé 
ation of flows on a steep slope begins in a region where the boundary laye! 
generated by the channel bed has reached the water surface, it is to be pre 
sumed that aeration depends primarily upon the turbulence intensity of the 
flow. Aeration appears to occur when the transverse velocities of turbule 
are sufficiently strong near the air-water interface to cause clumps of wa 
to break through the surface into the air and then fall back by gravity into 
flowing stream. This breaking away of clumps of water from the main st 
and the falling back into the stream with attendant splashing and breaking 
a heterogeneous spray of globules and droplets is associated with the inst 
flation of air into the stream. It occurs particularly with water flow at hi 
velocities on steep slopes. The air carried back is then distributed thro ul 
out the flow by turbulent transfer. Observations by means of instruments 
cate that two regions of self-aerated flow develop: (1) an upper region of 
heterogeneous clumps, globules, and droplets of water ejected from the fl 
ing liquid stream into the atmosphere at more or less arbitrary velocitie 


CE OPEN CHANNELS 1890-7 


1 (2) a lower region consisting of air bubbles distributed through the flow 
turbulent transport fluctuations, which can be described by some boundary- 
er equation. Between the two regions is a transition zone defined by a 
nsition depth which is a fluctuating surface necessarily at a statistical 

an elevation above the channel bottom. 


stribution of Air in Upper Region 
If it is assumed that the components normal to the bed of the turbulent 
ocity fluctuations are randomly distributed such that the mean value is 
ro, those in the outward direction will constitute one-half of a Gaussian 
tribution. Those in or near the surface that penetrate the surface will 
rry large clumps and smaller globules of water into the atmosphere against 
force of gravity a distance proportional to the square of the individual 
insverse velocities. The maximum height of the trajectories that they 
low will also be randomly distributed about zero. (The path will be a long, 
t trajectory because of the large, longitudinal component of the local ve- 
ity.) If the distances from the surface from which the particles are pro- 
ted are represented by one half of a Gaussian distribution, the frequency 
') of particles projected a distance y' above the transition depth may be ex- 
2ssed as 

(2)? 
iy} re seen : 

y har 


ere h is a measure of the mean distance the particles are projected above 
. The factor 2 is used to indicate that only particles being projected in the 
ward direction are being considered; that is, half of the complete Gaussian 


(4) 


co 
being taken so that , ‘i f(y') dy'=1. Furthermore, since it is presumed 


t the flow is uniform on the average, particles are equally likely to be pro- 
ted at all points of the transitional surface. Then the proportion of all 
‘ticles leaving a unit area of the transitional surface that reach or pass 
ough a corresponding area at elevation y' will be 


5 w sy ps 
alps 2 a wh 
Pye aa e dy (5) 


{the number of particles reaching or passing through this area in both di- 
tions during the sampling period will be 


Ny = N+ a (6) 


sre N,n is the total number of particles that leave or return to the transition 
face per unit area during the sampling period. The average number of 
ticles of water that reach or pass through any horizontal area per unit 

e may be assumed proportional to the average concentration of water parti- 
s. Therefore, using Eqs. (5) and (6), the air concentration at any distance 
\bove the transition level in terms of the concentration at the transition 

el Cp and the representative height h is 
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y'\2 
saa pera os 
I- Cy hr 4 y 


Eq. (7) is, of course, the cumulative Gaussian probability curve and is 4 ¥ 
cable to the aeration phenomenon only for y' 70. The gradient of the cont 
tration is from Eq. (7) 


ae 
dC 2c) -Gy) 
dy whim 


and when y' = 0, rs is a maximum whose value is 


aS 2(1-Cr) 
dy /max h./a 


and occurs where the concentration is C-p and the depth is dry. 

The constants Cp and h are properties of the flow related to the intens 
of the turbulence in the neighborhood of dp which, in turn, is probably de- 
pendent upon the velocity and roughness of the channel. 


Distribution of Air in Lower Region 

The air entrained at the transition surface by the turbulent fluctuations 
then distributed throughout the flow in the region below dp by the turbule 
in the stream resulting in a statistical equilibrium between the bouyancy ¢ 
the air and the concentration gradient. Such equilibrium is described by fl 
well-known equation 


‘dc 
dy 


-CWUte, =O 


where C is again the air concentration at a normal distance y from the bet 
Vp, the rising velocity of the air bubbles, is taken as being negative; and € 
is a mixing parameter for the air-bubble transfer. In aerated flows the Vv. 
of €) as a function of y is unknown, and a number of assumptions as to 
form of the function can be made. However, it may be presumed that the 
bubbles are transported in the same manner as momentum; and it appear 
reasonable as an approximation to assume €p to be proportional to €,,, th 
momentum mixing parameter, as is done in the case of transport of other 
substances such as sediment. In open channels in the nonaerated conditic 
where the shear is linearly distributed from zero at the surface to a maxX! 
mum at the bed and the velocity profile is considered to be logarithmic — 
throughout the depth of flow, it has been shown(9) that Em has a parabolic 
form being zero at the surface and at the bed with a maximum at the cent 
line and can be expressed in the form 


«5° Phyo (Ty 
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re VT ./p = V,, the shear velocity, and To is the boundary shear at the 

, P is the density of the fluid, 8 is a proportionality factor in Ep = Bem’ 

k is the Karman universal constant. In the case of transport in a closed 
nnel, however, the shear at the upper boundary is not zero and the maxi- 
m velocity occurs at a point between the upper and lower boundary, at 

ch point the shear is also zero. This requires that €m also be zero at this 
at and that the distribution of €,, consists of two parabolas having zero at 
bed, the point where the maximum velocity occurs, and at the upper 
ndary and maximum values at intermediate points above and below the 

at of maximum velocity. It has been shown,(10) however, that suspended 
iter is actually transported transversely across the level of zero mixing 
fficient, and it has been suggested that the mixing coefficient be considered 
stant in this region. Since it is likely that the vertical transport of sus- 
ded matter depends only upon the transverse velocity fluctuations in the 
tical direction, such transport may well occur even if the correlation be- 
en the vertical and longitudinal fluctuations approaches zero. 
Measurements of velocity distribution in aerated flows show(4,5) that the 
ocity tends to decrease as the upper limit of the flow is approached and 

t the maximum occurs well down in the flow proper. This implies that a 
ar is developed in the transition zone. A plausible explanation for such 
hear is the change in momentum engendered by the return of water clumps 
droplets to the flow at the end of their trajectory through the atmosphere 
ve dp. In the course of that trajectory, they suffer a loss of velocity due 
itmospheric drag on the mass. The shear developed at dp is analogous to 
ipper boundary and a distribution of the type found for closed channels 

tht be expected. 

The complex distribution found for closed channels may be approximated 

4 parabolic distribution over the entire region below dy, with the maximum 
ie incorporated in the factor 8B. This approximation is equivalent to as- 
ting a distribution for €p as given by Eq. (11). If this value for €p, given 
tq. (11), is substituted, Eq. (10) becomes 


Se Ass ek la le 
CV, = Bk 7p ( dy \y dy | (13) 


ch upon integration, assuming that Vj is independent of y, is 


« y 
evs C am (14) 


vyhich 


Vb Vp 
cai Ei hee 15) 
BkYtc/p BkVx ( 
dy 
Cy is a constant whose value is the concentration at y = 2 


| 
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Application to Experimental Data 


The use of Eqs. (7) and (14) as a description of the air-concentration 
distribution characteristics depends upon the degree to which they fit the 
perimental data. There are two factors which have a bearing on this: th 
form of the equation and the magnitude of the constants. Whether or no 
form of the equations is reasonably correct depends upon the extent that 
assumptions made in their development approximate the actual mechani, 
The constants involved depend upon the intensity of the turbulence and a 
present cannot be determined analytically but must be evaluated empiric; 
from systematic experimentation, 

Eq. (9) permits the determination of the transition depth dy and the co 
tration Cy at that depth. From a plot of the concentrations in terms oft 
normal distance from the bed, the point where the gradient i. is a maxin 
was graphically located, as shown in Fig. 7(a). The corresponding value 


dy and Cp are noted. Using this value of C-y, the ratio the ) was plot 


as a function of y on so-called probability paper, upon which a cumulatiy 
Gaussian distribution plots as a straight line. The value of h, or rathert 
value of o, the standard deviation of the distribution, wherein h -oy?2, is 
tained from the slope of this straight line. Alternately, h can be determi 
directly from Eq. (9) by evaluating the gradient at the point where = is 
maximum and inserting the value of Cp obtained for that point. Fig. 7(b) 
shows a plot of the concentrations in the upper region of the flow plotted 
this manner with the standard deviation indicated. It is apparent from tl 
plot that the concentrations possess a cumulative Gaussian distribution. 

In the region of flow below dy, Eq. (14) indicated that the concentratic 


is a function of eee which, when plotted logarithmically, results in a 


aT 
straight line whose slope is equal to the exponent z and Cj is the value 0 
at roms = 1. The concentration data for the experiment have been plot 


this manner in Fig. 7 to show the degree of agreement with the form of I 
(14) and to indicate the evaluation of Cy and z. Inasmuch as C approache 
infinity as y approaches dy, the experimental data must depart from the 


for the larger values of reves , but in this case the data agree with the 
up to 7s = 0.9. 


T 

The straight lines drawn through the data in Figs. 7(b) and 7(c) have | 
transposed to the original plot in Fig. 7(a) to show the characteristics ol 
curves on the arithmetic scales. Eqs. (7) and (14) have also been plotter 
along with the data for a number of typical experiments using various ¢ i 
charges and slopes shown in Figs. 8 and 9 for purposes of comparison. 
cause of the character of Eq. (14) in the neighborhood of y = dp, the two 
curves are connected by a transition to pass from one curve to the othe! 


Relation of Parameters to Flow Conditions 


Based upon the above analysis of the aeration phenomenon, a number 
parameters were developed as being descriptive of both the magnitude 4 
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tribution of the air. These parameters varied with the flow conditions 

h as the channel slope and discharge, which in turn, along with the rough- 
s characteristics of the channel, govern intensity and scale of the generat- 
turbulence. The flow turbulence is created initially at the bed by the wakes 
| eddies formed by the flow over the roughness elements; the turbulent 

lies are then diffused upward into the flow stream. The intensity of turbu- 
ce at the transitional surface which appears to be paramount in the aeration 
cess should then depend upon both the initial generation and the depth of 
mixture to the transitional surface. The intensity of the turbulence is re- 
2d to the boundary shear T, on the bed, which can be expressed by a so- 
led shear velocity Veep = VTo/p. Because of the shear that is presumed to 
st at the transitional surface, it is not possible in the present experiments 
lifferentiate the bed shear from the total shear. As a first approximation, 
total shear or a shear velocity based upon the total boundary shear may 
used as a measure of turbulence intensity. To test this hypothesis, the 

an concentration was plotted in terms of V,/ay2 3 where V,, is defined as 


Vy = /gd_sina (16) 


| dp is the transitional depth previously defined. The expression V,/dp2/ 3 
=mpirical and was chosen as that form which best correlated the data. 

'S plot is shown in Fig. 10. The measured values of the mean concentration 
_ reasonably well on a single curve with a tendency for the points corre- 
mding to the smaller discharges to depart from the mean curve. When the 
an concentration in the area below dy, Cr, is considered as being the air 
ch is insufflated and distributed through the flow, the correlation is even 
ter, as shown in Fig. 11. For the higher discharges for each slope the 

nts all fall on a smooth curve, while there is a clear departure from the 
ve for the low discharges on the corresponding slopes. It appears that as 
discharge decreases below a certain value for a given slope, an additional 
tor becomes important in reducing the concentration. The nature of this 
tor is not readily apparent, but it should be noted that a decrease in dis- 
rge corresponds to a decrease in dy, and it is quite possible that an insta- 
ty in flow pattern arises when the depth becomes small. Considering the 

Z as a whole, it is rather significant that when plotted in this manner, mean 
concentrations ranging from about 0.05 to 0.75 corresponding to discharges 
ging from about 4 to 15 cfs on slopes with ranges from 7.5 to 75 degrees 
correlated along a single curve. 

The parameter V. / Ap? 3 includes the depth dy to the transitional surface. 
s depth is a function of the air-concentration distribution and is difficult 
letermine explicitly. It is, however, related to the discharge and slope and 
‘eneral to the channel roughness. In these experiments the channel rough- 
s is a constant so that for these data, dy depends only on the slope and 
charge. In terms of these measured characteristics, V,/ay2 3 is related 
/q1/5 where S = sina and a is the angle of inclination of the bed and q is 
unit discharge of water. Consequently, the air concentration that may be 
ected in a high-velocity open channel can be related to the flow charac- 
stics. This relationship is shown in Fig. 12, where the mean concentration 
neasured is plotted as a function of S/qi/ 5 for the particular channel in 

ch the experiments were performed. For channels of different roughness 
racteristics, the curve will probably be similar but shifted an amount to 
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correspond to the different roughness. A similar plot is given in Fig. 13 
showing the mean concentration between the bed and the transitional surf 
where the depth is dy. Here again, as in Fig. 11, correlation between Cr 
and S/q1/5 appears to be better than that between C and S/q1 
The results shown in Fig. 13 provide a means of estimating the mean 
concentration that may be expected in an equally rough, wide open channel 
known unit discharge and slope. The depth of flow or the magnitude of th 
various previously described depth parameters, which are important ind 
sign, is a function of the air concentration. 


Effect of Entrained Air on the Flow 


It was suggested above that the air insufflated into the flow is distribu 
throughout the flow by turbulent transport. The entrained air has the eff 
of changing the flow from one of water alone to the flow of a mixture of a 
and water. Several characteristic depths were defined, such as d,, to rep 
sent the uppermost extent of the water particles projected from the flow, 
to represent the depth of flow of the water alone, and dp to represent af 
equilibrium depth brought about by the presence of air in the flow. Fore 
sign purposes, it is of considerable importance to know how these param 
and the actual flow velocity are related to the bulk flow properties or whi 
differences might be expected in these parameters in relation to a corre- 
sponding nonaerated flow for which the depth may be computed by means 
well-established flow formulas such as the Manning or Chezy formula. 

In order to investigate such a relationship, experiments were perform 
in the experimental channel to establish the resistance of the channel to! 
aerated flow, that is, flow of water alone. The slope was decreased so th 
for comparable discharges the flow would be nonaerated. These slopes 
ranged from approximately 1 degree to a maximum of about 6.5 degrees, 
above which aeration was observed to occur. The discharges were in the 
same range as those used in the aeration experiments. Normal flow was 
tablished for each combination of slope and discharge, and the mean dept 
was determined by averaging the depths measured at various points alor 
centerline. 

The relationship between the depth of flow for the nonaerated conditic 
the discharge and slope is shown in Fig. 14, where the measured depth d, 
plotted as a function of q/S 1/2 on logarithmic scales. The equation of th 
line drawn through the points is 


qs" 90sa reine 


For the range of the variables encountered in the nonaerated flow, it app 
that the Chezy formula best describes the data and the coefficient 90.5 ¢! 
terizes this particular roughness. The mean velocity is then 


Vm= 90.5dq “(sina)”? 


If it were assumed that the same type of flow equation applied to ae: a 
flow, the mean velocity of the air-water mixture flowing at a depth dip is 


CE OPEN CHANNELS 1890-13 


Marie deea(sin a) (19) 


ere V is the mean velocity of the air-water mixture and C is the Chezy 
efficient. Assuming that the velocity of the air and water components are 
ial, then Eq. (19) can be written in terms of the water discharge, as 


=, fe 

q = Cdd,’ “(sina)”? (20) 
ere d is the mean depth computed by Eq. (1) and represents the cross- 
stional area of the water. Eq. (20) was tested by plotting the experimental 
-a for aerated flow in the form ddpt 2 as a function of q/(sin o)1/2 in Fig. 
, Superimposed upon the plot is the line for nonaerated flow in the experi- 
tal channel taken from Fig. 14. The data cluster about the straight line 
*nonaerated flow, indicating that an equation similar to Eq. (20) will corre- 
e the data, and, in addition, when the parameters d and dy are used, the 
ezy resistance coefficient is essentially the same as that for nonaerated 
w. Fig. 15 shows that at least for the conditions of these experiments, air- 
‘rained flow, when considered to be a flow of a mixture, follows the same 
vs as for nonaerated flow when depth parameters related to the properties 
the mixture are used to define the linear flow dimensions. The agreement 
the data for aerated and nonaerated flow shown in Fig. 15 substantiate 
mewhat the deductions made in regard to the air-entrainment process. 

An evaluation of d, dy, and dy for specific flow conditions is not possible 
m Fig. 15. A clearer idea of the changes in these parameters will be ob- 
ned if they are compared with the depth d,, computed in the usual manner, 
ng Eq. (17) which includes for this channel the Chezy coefficient of 90.5. 
S the depth at which the water would flow in this channel at the given dis- 
irge and slope assuming no air entrainment. The ratios dy/dm, dp/dm, 
1 d/d,, as functions of mean air concentration are plotted in Fig. 16(a), 
b), and 16(c). In all cases, of course, the ratios approach unity as the 
an_concentration approaches zero. That is, for water flow only, j 
=d=dy=dm. As the mean air concentration increases, the ratios depart 
m unity in such a way that as C approaches, 1.0, d becomes very small, 
ile dy and dy necessarily increase. 
The upper limit of the spray d, increases very considerably with increased 
concentration which in turn means an increase in the turbulence intensity 
sroject water particles farther into the atmosphere. For the highest 
centrations measured, the value of d, is nearly four times as large as the 
aerated depth for a corresponding discharge and slope. The transition 
th, as might be expected, also increases relative to the computed depth 
ce increasing concentration bulks the flow and makes the effective depth 
ater. On the other hand, the mean depth decreases with increasing concen- 
tion. It appears from Fig. 16 that concentration has relatively little effect 
dp and d until the air concentrations are greater than about 50 per cent. 
>the slopes where the air concentrations were the greatest (being about 85 
cent), d was reduced to about one half of dm while dy was nearly twice 
arge as dy. Fs ‘ 

he decrease in d with increase in mean air concentration implies a corre- 
ding increase in the mean water velocity above that computed on the basis 
m Since d is the cross-sectional area of the water flow. The trend of the 
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mean water velocity in terms of the nonaerated velocity in this channel 
shown in Fig. 17, Here again, the velocity of aerated flow is about the sa 
as that of a corresponding nonaerated flow until the mean concentration is” 
the neighborhood of 50 per cent. It then increases rather rapidly. 


CONCLUSIONS 


1. Self-aerated, high-velocity, open-channel flow appears to consist of 
two distinct phases: an upper region consisting primarily of independent 
droplets and larger agglomerations of water that move independently of the 
stream proper, and a lower region in which discrete air bubbles are sus 
ed in a turbulent stream and are distributed by the mechanism of turbule 

2. In the upper region, the distribution of water-air agglomerate and 
lets can be adequately described by the Gaussian cumulative probability 
equation, the constants of which are properties of the flow conditions. In 
lower region, the distribution of the air agrees closely with an equation 01 
turbulent mixing based upon an approximation for the distribution of the m 
ing parameter. 

3. The magnitude of the entrained air concentrations apparently depend 
upon the intensity of turbulent fluctuations generated at the rigid bed andl 
depth of flow. The mean concentration could be correlated with the paran 
V,/dy2/3 or in terms of the flow characteristics; the mean concentration 
is a function of S/ql/ 5, where § is the sine function of the slope angle, a 
is the unit discharge. 

4, The maximum height of water (or spray) above the rigid bed relatii 
a corresponding nonaerated flow increases rapidly with mean air concen- 
tration in the flow. The transition depth is greater than the depth of a cc 
sponding nonaerated flow because of the bulking effect of the entrained ai 

5. The mean velocity of an air-entrained flow is greater than that of ; 
corresponding nonaerated flow by an amount that increases with the air 
concentration and corresponds to a decrease in the mean depth. ; 

6. When the effective depth of aerated flow is taken to be the transitie 
depth and the mean depth is the cross-sectional area per unit width of 
the flow is described by a formula of the Chezy type involving approxime 
the same roughness coefficient as applies for nonaerated flow in the sam 
channel. 
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LIST OF SYMBOLS 


( - Concentration, ratio of volume of air to volume of air plus water. 
a - Mean concentration in vertical section. 

C,;  - Concentration at point y = dy /2. 

Cy - Mean air concentration in region between lower boundary and dr. 
Cy - Concentration at transition depth y = dp. 

d - A mean depth of flow defined by Eq. (1). 

d, | - Upper limit of flow or value of y where Cc = 0.99. 

dp - Transition depth between upper and lower regions of flow. 

d,, ~ Mean depth in nonaerated flow. 

& 


- Acceleration due to gravity. 
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h - A mean height to which water particles are projected above « 
- Von Karman universal constant for velocity distribution. 


k 
Ny: - Number of particles reaching or passing through a horizo 
a distance y' above transition depth dr. 


Ny - Number of particles that leave or return to an area at the t 


depth. . 

Py: - Proportion of all particles leaving area at transition depth 
reach or pass through horizontal area at y' above transition 

Q - Total-water discharge. 

q - Water discharge per unit width of channel. 

V - Local velocity. 

Vv - Mean velocity in vertical section = q/d(d). 

V,  ~ Shear velocity = V g(dy) sina . 

Vp  - Rising velocity of air bubbles. 

y - Normal distance from channel bottom. 

y' - Outward normal distance above transition depth. 

Z - Exponent in air-concentration equation. 

a - Angle of inclination of channel. 

8B - Constant. 

Sb - Mixing coefficient for air bubbles in turbulent flow.. 

Em 7 Mixing coefficient for momentum transfer. 


- Kinematic viscosity of water. 
- Density of water. 
- Standard deviation of air distribution above dy. 


- Local shear force. 


oly A eS 


- Boundary shear force. 
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Fig. 6 - Air Concentrations as Functions of Slope and Discharge 
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Fig. 10 - Mean Concentration for Entire Flow as a Function of v/a, 
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Fig. 12 - Mean Concentration for Entire Flow as a Function of s/q/ 5 
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Fig. 13 - Mean Concentration in Region Below Transition Depth as a Function 
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Fig. 14 - Relationship of Depth to Discharge and Slope for Nonaerated 


Flow in Experimental Channel 
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Fig. 15 - Relationship of Depth to Discharge and Slope for Aerated Flow in 
Experimental Channel 
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Fig. 16 - The Influence of Air Concentration on Depths of Aerated Flow Relative 
to Corresponding Nonaerated Flow 
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December, 1958 


(Table Il - Summary of Experimental Measurements—Rough Channel) 


ra aa te Mere 2.2 021209 
201 Ook 3.2 0.162 
202 0 he2 0.1937 
203 (0) 562 0.220) 
20 ) 6.4 0.2515 
205 Ry) 7e2 002671 
206 ) A.2 -02291:3 
207 04 9.6  0.32h7 
910 15° 21.0 2.2 0.1010 
211 ) Sea OE 28O 
212 0 e200 5b 
213 0 5.2 0.1701 
21h ) 6.4 0.1937 
215 re) 7.2 022087 
216 re) 8.2 0.2275 
217 R) 9.6 0.2591 
220 22.5° 19.6 2.2 0.0886 
221 0.8 Jed 068230 
222 0 he2 0.1358 
223 ape pea Aare 
22h (e) 6.4 0.1763 
225 0.1 732 068929 
226 0.1 8.2 0.2226 
227 0.68 9.6 0.2329 
228 1.2 12.8 - 0.2855 
230 30.0° 21.0 2.2 0.0760 
231 2.2 3.2 0.0989 
232 n) Let 0.1325 
233 0 5.2 0.1307 
23 Tea 6.4 0.1530 
235 0.2 ee 016k 
236 0.2 8.2 0.1809 
237 1.0 9.6 0.2008 
238 0.2 dia »somee 
239 QO aT s0l2792 
2,0 als) 12.8 0.22) 


0.018 
0.0669 
0.0860 
0.1043 
0.127) 
0.1400 
0.1605 
0.1865 


0.0322 
0.0168 
9.9690 
9.0719 
0.0866 
0.0968 
9.1097 


0.0273 
0.0398 
9.0517 
0.0640 
9.0768 
0.0877 
0.1009 
0.1179 
9.1554 


0.022h 
0.0332 
0.0419 
0.0522 
0.0655 
0.0736 
0.0648 
0.0989 
0.1121 
0.12h1 
0.1271 


Table II - Summary of Experimental Measurements—Rough Channel 


ASCE 


OPEN CHANNELS 


psteer 


er eee ee ae” eee SS eee ae he ae 


| Ay?? Tre egaler 


0.45 
0.50 
0.50 
0.50 
0.50 
0.54 
0.48 
0.50 


0.008 
0.007 
0.905 
0.007 
0.005 
0.006 
0.00) 
0.006 


0.031 
0.047 
0.060 
0.067 
0.083 
0.072 
0.967 
0.058 


0.107 
0.115 
0.123 
0.116 
0.140 
0.123 
0.107 
0.099 
0.963 


0.222 
0.20 
0.245 
0.242 
0.213 
0.227 
0.200 
0.186 
0.185 
0.157 


0.183 


12.13 
12.97 
14.45 
15-75 


16.95" 


17.97 
18.55 
19.71 


14.50 
16.53 
18.50 
20.40 
22.0 
23.0 
24.0 
2h.7 


16.55 
18.85 
20.6 
221 
24.1 
2h.9 
25.8 
2765 
Coe) 


19.30 
21.7 
24.9 
26.5 
2769 
292 
30.2 
31.9 
3309 
33.8 
3502 


0.710 
0.809 
0.909 
0.969 
1.039 
1.07 
1.19 
1.178 


0.923 
1.050 
1.144 
1.206 
1.292 
1.340 
1.396 
1.49) 


1.081 
1.235 
1.338 
1.430 
1.530 
1.600 
1.676 
1.768 
1.911 


2.92 565 
2067 4.52 
2.68 5.91 
2.62 7.16 
2256 9205 
2.54 10.23 
eG 

2.46 13.40 
yl 5.58 
4.04 4.53 
3.92 6.03 
3.86 7034 
SLT 9.14 
3273 10.23 
3.63 11.59 
3.60 13.55 
5019 e771 
4.96 4.98 
By 6.21 
Ws72 7258 
4.62 9.54 
4.56 10.75 
eld Lose 
441 14.77 
4.30 19.10 
6.00 6.33 
Boe 52 
5655 7.18 
5 l3 8.79 
5.31 11.04 
5025 12.23 
5.16 13.83 
5.09 16.30 
5.02 18.36 
5.00 18.98 
4.96 20.80 
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1890-34 HY 7 . December, 1958 


[see [es [ni [aL ee 


20} 37.5° 19.6 0.9605 0.198 0.985 0.9176 
21 


16.8 p is 0.0783 0.245 0.120 0.0271 2 
22 16.8 2 0.0923 9.291 0.148 0.0355 3.59 
2h)3 16.8 52 0.1073 0.326 0.175 O.0b49 hekS 
ahh 5.0 64 0.1317 0.361 0.195 0.0582 5.46 
2hs 10.6 ee. 0.149, 0.388 0.212 0.0688 6.15 
2h6 93 Bea 0.1638 0.413 0.227 0.9780 7.01 
27 10.5 9.6 0.1855 0.433 03250 0.0927 8.21 
28 10.3 SHEP: 0.206). 0472 0.272 0.1075 9657 
2h9 10.3 12.8 002275 O«51h 0.300 94126 10.93 
250 10.3 15.9 0.2521 0.505 0.325 0.1415 12.81 
2503 5° 18.5 2.2 0.0530 0.188 0.090 0.0159 1.7h 
251 12.9 3.2 0.0625 0.2h0 0.125 0.0221 2.53 
252 12.8 he2 0.9791 0.303 0.151 0.0307 3.33 
253 12.8 5.2 0.0921 0.342 0.172 0.0382 4.13 
25h 12.9 6.4 0.1065 0.378 90.205 0.0482 5.07 
255 12.0 Tee 0.1182 0.400 0.220 0.0555 5.71 
256 TIA 8.2 0.1361 0.422 0.20 0.0667 6.50 
257 11.9 9.6 0.1600 0.452 0.260 0.0815 7.61 
258 JES 2IS2 0.1782 0.462 0.260 0.0910 8.87 
259 18.5 12.8 002029 O06h79 02275 0.1063 10.2) 
260 19.8 15.0 0.2397 0.469 0.278 0.1280 11.90 
261 60° 19.8 2.2 0.0397 0.190 0.085 0.0116 1.58 
262 14.5 3.2 0.0540 0.260 0.126 0.0195 2.29 
263 18.0 4.2 0.0648 0.290 0.160 0.0259 3.01 
26), 18.6 52 9.0759 0.331 0.190 0.0330 3.73 
265 18.6 604 0.0887 0.389 0.220 0.0416 4.59 
266 19.8 Teo 0.1003. 0.423) (0.20) 0.01920 Soak 
267 19.2 8.2 0.1163 0.476 0.260 0.0593 5.88 
268 14.5 9.6 0.1299 0.491 0.280 0.0687 6.88 
270 7e° S298 Dee 0.03@2 0.212 0.100 0.0121 1.h9 
271 22.0 3.2 0.0429 0.273 0.140 0.0161 2.17 
272 22.0 4.2 0.0618 0.333 0.200 0.0275 2.85 
273 eles 52 0.0764  O.440 0.220 0.0358 3.53 
27h 21.3 6.4 0.087) 0.501 0.250 0.0457 h.35 
275 20.9 Te2 0.0936 0.465 0.275 0.0490 .88 
276 20.9 8.2 0.1016 0.500 0.270 0.0328 5.56 
277 20.9 9-6 0.112 0.520 0.290 0.0606 6.51 
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‘ 0.63 0.431 
3 «60.67 = (050 
> 0.68 0.42 
0.69 0.72 
>? 0.68 0.402 
5 60 0662020392 
L 0.63 0.363 
mm 0.64 0.317 
Mm 0.62 0.310 
) 0.6 0.292 
mm 0.61 0.22 
0.72 0.510 
5 0.76 0.565 
) 0.75 0.550 
> «600673 = (06555 
Me O72 02555 
me 0.73 «860.525 
> ©6072): 0 90 
60.700 | (0.53 
2 6 0665S «(0398 
2 0.6 0.326 
0.60 0.217 
im 0.78 0.613 
) 0.80 0.630 
m 0.80 0.545 
2 60682) (065 
» O81 0.657 
1 0.80 0.640 
3 0.78 0.610 
»- 0.77 0.590 
? 0.83 0.690 
me 0.87 «6 0750 
m™ 0.85 0.730 
me 0.04 0.700 
mz O0.82 0.715 
m0.83 0.701 
im 0.81 0.687 
1 0.80 0.679 


OPEN CHANNELS 


6.66 Tehh9 
6.33 11.10 
6.09 15.23 
5.94 19.13 
5.82 15.h0 
5.73 19.60 
5.67 21.0 
5259 2.8 
5.51 26.1 
52 Sal 
5.38 35-7 
7.10 8.80 
6.73 12.60 
6.52 16.31 
6.38 19.75 
6.21 2563 
6.15 2667 
206 29.0 
5.98 31.0 
5.98 38.6 
5.91 41.0 
5.e6 4.0 
7293 11.50 
748 15.61 
7215 24.1 
6.95 30.8 
6.84 37.6 
6.68 42.1 
6.59 44.0 
6.50 4he2 
8.18 14.32 
7015 26.7 
7230 3.6 
7.18 3767 
7-03 46.1 
6.93 52.3 
6.93 
6.8 
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0.164 0.08 
0.152 0.056 
0.142 0.061 
0.10 0.068 
0.177 0.073 
0.158 0.078 
0.192 0.081 
0.250 0.080 
0.237 0.090 
0.290 0.089 
0.318 0.080 


0.140 0.046 
0.113 0.051 
0.113 0.068 
0.087 0.078 
0.092 0.079 
0.113 0.082 
0.122 0.080 
0.143 0.090 
0.151 0.088 
0.234 0.095 
0.346 0.083 


0.107 0.053 
0.102 0.061 
0.091 0.061 
0.098 0.073 
0.0806 0.090 
0.086 0.095 
0.087 0.098 
0.093 0.102 


0.093 0.060 
0.089 0.067 
0.063 0.084 
0.067 0.100 
0.050 0.125 
0.0579 0.09 
0.0506 0.118 
0.05075 0.106 
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The technical papers published in the past year are identified by number below. chnic 
sponsorship is indicated by an abbreviation at the end of each Paper Number, the symbols referring t€ 
Transport (AT), City Planning (CP), Construction (CO), Engineering Mechanics (EM), Highway ( 
draulics (HY), Irrigation and Drainage (IR), Pipeline (PL), Power (PO), Sanitary Engineering (8A) 
Mechanics and Foundations (SM), Structural (ST), Surveying and Mapping (SU), and Waterways and | 
(WW), divisions, Papers sponsored by the Department of Conditions of Practice are identified by t € 8 
(PP). For titles and order coupons, refer to the appropriate issue of “Civil Engineering.” B zinnin 
Volume 82 (January 1956) papers were published in Journals of the various Technical Divisions. To 
papers inthe Journals, the symbols after the paper number are followed by a numeral designating ti 
of a particular Journal in which the paper appeared. For example, Paper 1859 ie identified as 1858) 
which indicates that the paper is contained in the seventh issue of the Journal of the Hydraulie 
during 1958. 


VOLUME 83 (1957) 


DECEMBER: 1449(HY6), 1450(HY6), 1451(HY6), 1452(HY6), 1453(HY6), 1454(HY6), 1455(HY6), 14560) 
1457(PO8), 1458(PO8), 1459(PO6), 1460(PO6)°, 1461(SA6), 1462(SA5), 1463(SA6), 1464(SA6), 1466 
1466(SA6)©, 1467(AT2), 1468(AT2), 1469(AT2), 1470(AT2), 1471(AT2), 1472(AT2), 1473(AT2), 1476 Ae 
1475(AT2), 1476(AT2), 1477(AT2), 1478(AT2), 1479(AT2), 1480(AT2), 1481(AT2), 1482(AT2), 14885) 
1484(AT2), 1485(AT2)°, 1486(BD2), 1487(BD2), 1488(P06), 1489(PO6), 1490(BD2), 1491(BD2), 1492 
1493(BD2), one 


in 


VOLUME 84 (1958) o 


JANUARY: 1494(EM1), 1495(EM1), 1496(EM1), 1497(TR1), 1498(IR1), 1499(1R1), 1500(1R1), 1501(URE), 
(IR1), 1503(IR1), 1504(IR1), 1505(IR1), 1506(IR1), 1507(1R1), 1508(ST1), 1509(ST1), 1510(S8T1), 1921 
1512(ST1), 1513(WW1), 1514(WW1), 1515(WW1), 1516(WW1), 1517(WWI1), 1518(WWi), 1519(8TT) 
(EM1)°, 1521(1R1)°, 1522(ST1)°, 1523(WW1)°, 1524(HW1), 1525(HW1), 1526(HW1)°, 1527(Wi). 


FEBRUARY: 1528(HY1), 1529(PO1), 1530(HY1), 1531(HY¥1), 1592(H¥1), 1539(SA1), 1534(SA1), 1595 
1536(SM1), 1537(SM1), 1538(POL)°, 1539(SA1), 1540(SA1), 1541(SA1), 1542(SA1), 1543(SA1), 15444 
1545(SM1), 1546(SM1), 1547(SM1), 1548(SM1), 1549(SM1), 1550(SM1), 1551(SM1), 1552(SM1), ; 
1554(PO1), 1555(PO1), 1556(PO1), 1557(SA1)°, 1558(HY¥1)°, 1559(SM1)°, Mee, 


MARCH: 1560(ST2), 1561(ST2), 1562(ST2), 1563(ST2), 1564(ST2), 1565(ST2), 1666(ST2), 1567(ST2) 
(WW2), 1569(WW2), 1570(WW2), 1571(WW2), 1572(WW2), 1573(WW2), 1574(PL1), 1575(PL1), 76 
1577(PL1), 1578(PLi)°, 1579(ww2)°. i: 


APRIL: 1580(EM2), 1581(EM2), 1582(HY¥2), 1583(HY2), 1584(HY¥2), 1585(HY2), 1586(HY2), 1587(H¥2), 
(HY2), 1589(TR2), 1590(IR2), 1591(IR2), 1592(SA2), 1593(SU1), 1594(SU1), 1595(SU1), 1596(EM2), 1507 
1598(PO2), 1599(PO2), 1600(PO2), 1601(PO2), 1602(P02), 1603(HY¥2), 1604(EM2), 1605 
1607(SA2), 1608(SA2), 1609(SA2), 1619(SA2), 1611(SA2), 1612(SA2), 1613(SA2), 1614(SA2)°, 1615¢ 
(HY2)€, 1617(SU1), 1618(PO2)°, 1619(EM2)°, 1620(CP1). 


MAY: 1621(HW2), 1622(HW2), 1623(HW2), 1624(HW2), 1625(HW2), 1626(HW2), 1627(HW2), 1626 
(ST3), 1630(ST3), 1631(ST3), 1632(ST3), 1633(ST3), 1634(ST3), 1635(ST3), 1636(ST3), 1637(ST3) 
1639(WW3), 1640(WW3), 1641(WWS), 1642(WW3), 1643(WW3), 1644(WWS), 1645(SM2), 1646( 
(SM2), 1648(SM2), 1649(SM2), 1650(SM2), 1651(HW2), 1652(HW2), 1653(Wws3)°, 1654(SM2), 1656 
1656(ST3)°, 1657(SM2)°, ie 

JUNE: 1658(AT1), 1659(AT1), 1660(HY3), 1661(HY3), 1662(HY3), 1663(HY3), 1664(HY¥3), 1665(SA3 
(PL2), 1667(PL2), 1668(PL2), 1669(AT1), 1670(PO3), 1671(PO3), 1672(POS), 1673(PL2), 1674(PL2 
(P03), 1676(PO3), 1677(SA3), 1678(SA3), 1679(SA3), 1680(SA3), 1681(SA3), 1682(SA3), 1683(POS), & 
(HY3), 1685(SA3), 1686(SA3), 1687(PO3), 1688(SA3)°, 1689(P03)°, 1690(HY3)°, 1691(PL2)°. 


JULY: 1692(EM3), 1693(EM3), 1694(ST4), 1695(ST4), 1696(ST4), 1697(SU2), 1698(SU2), 1699(SU2), 1 
1701(SA4), 1702(SA4), 1703(SA4), 1704(SA4), 1705(SA4), 1706(EM3), 1707(ST4), 1708(ST4), 1709(STS), 
(ST4), 1711(ST4), 1712(ST4), 1713(SU2), 1714(SA4), 1715(SA4), 1716(SU2), 1717(SA4), 1718(EM3) 
(EM3), 1720(SU2), 1721(ST4)°, 1722(ST4), 1723(ST4), 1724(EM3)°, at 

AUGUST: 1725(H¥4), 1726(HY4), 1727(SM3), 1728(SM3), 1729(SM3), 1730(SM3), 1731(SM3), 173218 
(P04), 1734(P04), 1735(PO4), 1736(PO4), 1737(PO4), 1738(PO4), 1739(PO4), 1740(PO4), 17416 
(P04), 1743(PO4), 1744(PO4), 1745(PO4), 1746(PO4), 1747(PO4), 1748(PO4), 1749(PO4), 


SEPTEMBER: 1750(1R3), 1751(1R3), 1752(1R3), 1753(QR3), 1754(R3), 1755(ST5), 1756(STS), 275 
1758(ST5), 1759(ST5), 1760(ST5), 1761(ST5), 1762(STS), 1763(ST5), 1764(STS), 1765(WW4), F 
1767(WWw4), 1768(WWa4), 1769(WW4), 1770(Ww4), 1771(WW4), 1772(WWw4), 1773(Ww4), 17740 
(IR3), 1776(SA5), 1777(SA5), 1778(SA5), 1779(SA5), 1780(SA5), 1781(WW4), 1782(SA5), 1793/5 
(IR3)¢, 1785(WW4)°, 1786(SA5)°, 1787(STS)°, 1788(tR3), 1789(WW4). 


OCTOBER: 1790(EM4), 1791(EM4), 1792(EM4), 1793(©M4), 1794(EM4), 1795(HW3), 1796 (HWS), 
1798(HW3), 1799(HW3), 1800(HW3), 1801(HW3), 1802(HW3), 1803(HW3), 1804(HW3), 1805 
{HY5), 1807(HY5), 1808(HY5), 1809(HY5), 1810(HY5), 1811(HY5), 1812(SM4), 1813(SM4), 1814¢ 
{ST6), 1816(ST6), 1817(ST6), 1818(ST6), 1819(ST6), 1820(STS), 1821(ST6), 1822(@M4), 18230) 
(SM4), 1825 (SM4), 1826(SM4), 1827(ST6)°, 1828(SM4)°, 1829(HW3)°, 1830(PO5)°, 1631(BM4)°, 185 


NOVEMBER: 1833(HY6), 1834(HY6), 1895(SA6), 1836(ST7), 1837(ST7), 1838(ST7), 1839(ST7), 1B4l 
1841(ST7), 1842(SU3), 1843(SU3), 1844(SU3), 1845(SU3), 1846(SU3), 1847(SA6), 1848(SA6), 184 
Anh FAG) 1851(SA6), 1852(SA6), 1853(SA8), 1854(S7T7), 1855(A6)°, 1856(HY6)°, 1857(S7 
$U3)¢, 

DECEMBER: 1859(HY7), 1860(0R4), 1861(IR4), 1862(IR4), 1863(6M6), 1864(SM5 266 
(ST8), 1868(PP1), 1869(PP1), 1870(PP1), 1871(PP1), pili ivrg i Bovmasig yet err ' 
(WW5), 1877(CP2), 1878(ST8), 1879(ST8), 1880(HY7)°, 1881(SM5)°, 1882(ST8)° 1883(PPL)®, 1 84 
1885(CP2)°, 1886 (POS), 1887(PO6), 1888(PO6), 1889(PO6), 1890(HY7), 1891(PP1), a 


c, Discussion of several papers, grouped by divisions, ' ‘ 


